EFFECTS OF FRUSTRATED SURFACE IN HEISENBERG THIN FILMS 
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We study by extensive Monte Carlo (MC) simulations and analytical Green function (GF) method 
effects of frustrated surfaces on the properties of thin films made of stacked triangular layers of atoms 
bearing Heisenberg spins with an Ising-like interaction anisotropy. We suppose that the in-plane 
surface interaction J s can be antiferromagnetic or ferromagnetic while all other interactions are 
ferromagnetic. We show that the ground-state spin configuration is non linear when J s is lower 
than a critical value J^. The film surfaces are then frustrated. In the frustrated case, there are two 
phase transitions related to disorderings of surface and interior layers. There is a good agreement 
between MC and GF results. In addition, we show from MC histogram calculation that the value 
of the ratio of critical exponents j/v of the observed transitions is deviated from the values of two 
and three Ising universality classes. The origin of this deviation is discussed with general physical 
arguments. 
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I. INTRODUCTION 

This paper deals with the effect of the frustration in 
magnetic thin films. The frustration is known to cause 
spectacular effects in various bulk spin systems. Its ef- 
fects have been extensively studied during the last decade 
theoretically, experimentally and numerically. Frus- 
trated model systems serve not only as testing ground for 
theories and approximations, but also to compare with 
experiments^ 

On the other hand, surface physics and systems of 
nanoscales have been also enormously studied during the 
last twenty years. This is due in particular to applica- 
tions in magnetic recording, let alone fundamental the- 
oretical interests. Much is understood theoretically and 
experimentally in thin films where surfaces are 'clean' 
i.e. no impurities, no steps etc£&2£&l Less is known at 
least theoretically on complicated thin films with special 
surface conditions such as defects^, arrays of dots and 
magnetization reversal phenomenon. 10 ! 11 ' 12 ! 13 ' 14 ! 15 ' 16 

In this paper we study the frustration effect on prop- 
erties of thin films made of stacked triangular lattices. 
In-plane interaction of the surfaces is antiferromagnetic 
and that of interior layers is ferromagnetic. The film sur- 
faces are frustrated. 

The paper is organized as follows. Section II is de- 
voted to the description of our model. The ground state 
in the case of classical spins is determined as a function 
of the surface interaction. In section III, we consider the 
case of quantum spins and we apply the Green function 
technique to determine the layer magnetizations and the 
transition temperature as a function of the surface in- 
teraction. The classical ground state determined in sec- 
tion II is used here as starting (input)configuration for 



quantum spins. We are interested here in the effect of 
magnetic frustration on magnetic properties of thin films. 
A phase diagram is established showing interesting sur- 
face behaviors. Results from Monte Carlo simulations 
for classical spins are shown in section IV and compared 
to those obtained by the Green function method. Wc 
also calculate by Monte Carlo histogram technique the 
critical behavior of the phase transition observed here. 
Concluding remarks are given in Section V. 



II. MODEL 

It is known that many well-established theories failed 
to deal with frustrated spin systems i Among known 
striking effects due to frustration, let us mention the high 
ground-state (GS) degeneracy associated often with new 
symmetries which give rise sometimes to new kinds of 
phase transition. One of the systems which are most 
studied is the antiferromagnetic triangular lattice. Due 
to its geometry, the spins are frustrated under nearest- 
neighbor (NN) antiferromagnetic interaction. In the case 
of Heisenberg model, the frustration results in a 120° GS 
structure: the NN spins form a 120° angle alternately 
in the clockwise and counter-clockwise senses which are 
called left and right chiralities. 



A. Hamiltonian 

In this paper we consider a thin film made up by stack- 
ing N z planes of triangular lattice of L x L lattice sites. 
The Hamiltonian is given by 
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H = - 2J Ji,j$i ■ Sj - X! 

<*,3> 

where Si is the Heisenberg spin at the lattice site i,Yl{i j) 
indicates the sum over the NN spin pairs Si and Sj . The 
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last term, which will be taken to be very small, is needed 
to make the film with a finite thickness to have a phase 
transition at a finite temperature in the case where all 
exchange interactions Jij are ferromagnetic. This guar- 
antees the existence of a phase transition at finite temper- 
ature, since it is known that a strictly two-dimensional 
system with an isotropic non-Ising spin model (XY or 
Heisenberg model) does not have long-range ordering at 
finite temperature ^ 

Interaction between two NN surface spins is equal to 
J s . Interaction between layers and interaction between 
NN in interior layers are supposed to be ferromagnetic 
and all equal to J = 1 for simplicity. The two surfaces of 
the film are frustrated if J s is antiferromagnetic (J s < 0). 



B. Ground state 

In this paragraph, we suppose that the spins are classi- 
cal. The classical ground state (GS) can be easily deter- 
mined as shown below. Note that for antiferromagnetic 
systems, even for bulk materials, the quantum GS can- 
not be exactly determined. The classical Neel state is 
often used as starting configuration for quantum spins. 
We will follow the same line hereafter. 

For J s > (ferromagnetic interaction), the magnetic 
GS is ferromagnetic. However, when J s is negative the 
surface spins are frustrated. Therefore, there is a compe- 
tition between the non collinear surface ordering and the 
ferromagnetic ordering due to the ferromagnetic interac- 
tion from the spins of the beneath layer. 

We first determine the GS configuration for I = I s = 
0.1 by using the steepest descent method : starting from 
a random spin configuration, we calculate the magnetic 
local field at each site and align the spin of the site in 
its local field. In doing so for all spins and repeat until 
the convergence is reached, we obtain in general the GS 
configuration, without metastable states in the present 
model. The result shows that when J s is smaller than a 
critical value J£ the magnetic GS is obtained from the 
planar 120° spin structure, supposed to be in the XY 
plane, by pulling them out of the spin XY plane by an 
angle f3. The three spins on a triangle on the surface form 
thus an 'umbrella' with an angle a between them and an 
angle (3 between a surface spin and its beneath neighbor 
(see Fig. [I}. This non planar structure is due to the 
interaction of the spins on the beneath layer, just like an 
external applied field in the z direction. Of course, when 
J s is larger than one has the collinear ferromagnetic 
GS as expected: the frustration is not strong enough 
to resist the ferromagnetic interaction from the beneath 
layer. 

We show in Fig. [5]cos(a) and cos(/3) as functions of J s . 
The critical value J° is found between -0.18 and -0.19. 
This value can be calculated analytically by assuming 
the 'umbrella structure'. For GS analysis, it suffices to 
consider just a cell shown in FigJT] This is justified by 
the numerical determination discussed above. Further- 




FIG. 1: Non collinear surface spin configuration. Angles be- 
tween spins on layer 1 are all equal (noted a), while angles 
between vertical spins are f3. 
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FIG. 2: cos(a) (diamonds) and cos(/3) (crosses) as functions 
of J s - Critical value of is shown by the arrow. 

more, we consider as a single solution all configurations 
obtained from each other by any global spin rotation. 

Let us consider the full Hamiltonian . For simplicity, 
the interaction inside the surface layer is set equal J s 
(— 1 < J s < 1) and all others are set equal to J > 0. 
Also, we suppose that Zjj = I s for spins on the surfaces 
with the same sign as J s and all other Zjj are equal to 
I > for the inside spins including interaction between 
a surface spin and the spin on the beneath layer. 

The spins are numbered as in Fig. [1] Si, S 2 and S3 
are the spins in the surface layer (first layer), S[, S' 2 and 
S' 3 are the spins in the internal layer (second layer) . The 
Hamiltonian for the cell is written as 

H p = -6 [J, (Si • S 2 + S 2 • S 3 + S 3 • Si) 

+/ s (SISI + S%SI + SIS?) 
+ J (S[ • + • S 3 + S 3 ■ S'i) 

+1 {S[ Z S' 2 Z + S' 2 Z S' 3 Z + SgS'f)] 
- 2J(S 1 -S 1 + S 2 -S 2 + S 3 -S 3 ) 

-21 (SfS[ z + S' 2 Z S' 2 Z + S^S'i ) , (2) 

Let us decompose each spin into two components: an 
xy component, which is a vector, and a z component 

Si = (sj,Sf). Only surface spins have xy vector com- 
ponents. The angle between these xy components of NN 
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surface spins is which is chosen by (7^ is in fact the 
projection of a defined above on the xy plane) 



2-7T 47T 

71,2 = 0, 72,3 = -y ; 73,i = -g-- 



(3) 



The angles and /3- of the spin Si and S ■ with the z 
axis are by symmetry 



Pi = fa =03 = 13, 



The total energy of the cell with Si = = |, can 
be rewritten as 



9</s ■ 2 n 



HJs+Is) 2o 

2 C0S & 

(4) 



By a variational method, the minimum of the cell energy 
corresponds to 

—?-= (—J S + 9I S \ cos/3 sin/3 + ^(J + 1) sin/3 = 

(5) 

We have 



cos /3 : 



J + J 
9 J s + 6J S 



(6) 



For given values of I s and /, we see that the solu- 
tion (JSJ) exists for J s < J% where the critical value 
is determined by — 1 < cos/3 < 1. For / = — I s =0.1, 
Jg « —0.1889 J in excellent agreement with the numeri- 
cal result. 

The classical GS determined here will be used as input 
GS configuration for quantum spins considered in the 
next section. 



III. GREEN FUNCTION METHOD 

Let us consider the quantum spin case. For a given 
value of J s , we shall use the Green function method to 
calculate the layer magnetizations as functions of tem- 
perature. The details of the method in the case of non 
collinear spin configuration have been given in Ref. 18 . 
We briefly recall it here and show the application to the 
present model. 



A. Formalism 

We can rewrite the full Hamiltonian fTJ) in the local 
framework of the classical GS configuration as 



n 



<i,j> 



1 



(cos6 ij -i)(S+S+ + SrS7) 



+ i (cos ij + i) (stsi + srsf) 

+ \ sin% (5+ + Sr) SJ - \ sin%S* (S+ + S-) 



) <i,j> 



(7) 



where cos (9ij) is the angle between two NN spins deter- 
mined classically in the previous section. 

Following Tahir-Kheli and ter Haar^ we define two 
double-time Green functions by 



Gij(t,t') = «S+(i);S~(i')», 



(8) 
(9) 



The equations of motion for Gij(t, t') and Fij(t, t') read 
ijG itj (t,t') = ([St(t),Sr(t')])8(t-t>) 



(([H,St(t)];Sr(t'))), (10) 



z-F^M') - ([Sr(t),S-(t')])S(t-t') 



«[W,5r(*)];S-(f))), (ll) 



We neglect higher order correlations by using the 
Tyablikov decoupling scheme 2 ^ which is known to be 
valid for exchange termsj 2 ^ Then, we introduce the 
Fourier transforms 



Gij (t, t') — — 



dkxy 2ir 



j(t-t') 



Fi,j (*, t') 



9n,n' (^- , ? k-xy ) 6 



oo 

£ka. y -(Ri-Rj) 



(12) 



A J J dKv 2ir 

fn.n' (^5 ^-xy) & 



due 

oo 



('-*'). 



(13) 



where uj is the spin- wave frequency, \s. xy denotes the wave- 
vector parallel to xy planes, is the position of the spin 
at the site i, n and n' are respectively the index of the 
layers where the sites i and j belong to. The integral over 
k^j, is performed in the first Brillouin zone whose surface 
is A in the xy reciprocal plane. 

The Fourier transforms of the retarded Green functions 
satisfy a set of equations rewritten under a matrix form 



M(w)g= u, 



(14) 



where M (uj) is a square matrix (2N Z x 2N Z ), g and u 
are the column matrices which are defined as follows 



/ 9i,n> \ 

fl,n> 



9N z ,n' 



/ 2 <5f \ 





2(5^ )S N „ 




(15) 



4 



and 
M{uj) = 

where 
At 



( At B X 



-D7 



C 

-c 



N, 



c 

-c 



N, 



N, 



Bn z 



(16) 



-J n {S* n ) (Zi) (cos^ + 1) 



J n (S*) ZcOS0 n - J n ,n+1 (Sn+l) COS0 n>n+ i 
Jn,n—1 \£>n—l) cos ^n,ri-l — ZI n (S^) 



I, 



n,n+l 



B n = -J„(S*)(cosfl n -l)(Z7), 

Ct = \jn,n-l{S Z n ){cOsQ n , n - X ±l), 



J n ,n+1 (Sn) (C0s6» n , n+ i ± 1) 



(17) 
(18) 
(19) 
(20) 



in which. Z = 6 is the number of in-plane NN, Q n ,n±i 
the angle between two NN spins belonging to the layers 
n and nil, 9 n the angle between two in-plane NN in 
the layer n, and 
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2 cos (k x a) + 4 cos (k y a/2) cos (k y aV3/2\ jZ 



Here, for compactness we have used the following nota- 
tions: 

i) J n and /„ are the in-plane interactions. In the 
present model J n is equal to J s for the two surface layers 
and equal to J for the interior layers. All I n are set to 
be I. 

ii) Jn.n±i and I n n ±i are the interactions between a 
spin in the n th layer and its neighbor in the (n i l) th 
layer. Of course, J n , n -i = l n ,n-i=0 i£ n = 1, J n ,n+i = 
I n ,n+i=0 if n = N z . 

Solving det|M| = 0, we obtain the spin- wave spectrum 
lj of the present system. The solution for the Green func- 
tion g n<n is given by 



1 M L 

IMI ' 



(21) 



where |M| n is the determinant made by replacing the 
n-th column of |M| by u in (fTS"]) . Writing now 



(22) 



one sees that oj% (k Ky ) , i = 1, • • • , N z , are poles of the 
Green function g n ,n- O^xy) can be obtained by solving 
|M| = 0. In this be expressed as 



(23) 



where h n (uJi (k xy )) is 



h n (Wi (k xv )) 



|M|„ (Ui (Ky)) 



(24) 



Next, using the spectral theorem which relates the cor- 
relation function (S~ Sj~) to the Green functions^ one 
has 



(S7St) = lim 



1 

o A 



rfk 



2tt 



{gn,n> (w + ie) 



9n,n> ( w -* £ ))- 



,f)u _ 1 



(25) 



where e is an infinitesimal positive constant and (3 = 
1/ksT, ks being the Boltzmann constant. For spin S = 
1/2, the thermal average of the z component of the i — th 
spin belonging to the n — th layer is given by 



(SI 



= 2 



(26) 



In the following we shall use the case of spin one-half. 
Note that for the case of general S, the expression for 
(Sf) is more complicated since it involves higher quanti- 
ties such as ((S^) 2 ). 

Using the Green function presented above, we can 
calculate self-consistently various physical quantities as 
functions of temperature T. The first important quan- 
tity is the temperature dependence of the angle between 
each spin pair. This can be calculated in a self-consistent 
manner at any temperature by minimizing the free energy 
at each temperature to get the correct value of the angle 
as it has been done for a frustrated bulk spin systems^ 
In this paper, we limit ourselves to the self-consistent cal- 
culation of the layer magnetizations which allows us to 
establish the phase diagram as seen in the following. 

For numerical calculation, we used / = 0.1J with J = 
1. For positive J s , we take I s = 0.1 and for negative 
J s , we use I s = — 0.1. A size of 80 2 points in the first 
Brillouin zone is used for numerical integration. We start 
the self-consistent calculation from T = with a small 
step for temperature 5 x 10 -3 or 10 _1 (in units of J/ks)- 
The convergence precision has been fixed at the fourth 
figure of the values obtained for the layer magnetizations. 



B. Phase transition and phase diagram of the 
quantum case 

We first show an example where J s = —0.5 in Fig. [3] 
As seen, the surface-layer magnetization is much smaller 
than the second-layer one. In addition there is a strong 
spin contraction at T = for the surface layer. This is 
due to the antiferromagnetic nature of the in-plane sur- 
face interaction J s . One sees that the surface becomes 
disordered at a temperature T\ ~ 0.2557 while the sec- 
ond layer remains ordered up to Ti ~ 1.522. Therefore, 
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the system is partially disordered for temperatures be- 
tween T\ and T 2 . This result is very interesting because 
it confirms again the existence of the partial disorder in 
quantum spin systems observed earlier in bulk frustrated 
quantum spin systems i 18 i 23 Note that between T\ and 
T2, the ordering of the second layer acts as an external 
field on the first layer, inducing therefore a small value of 
its magnetization. A further evidence of the existence of 
the surface transition will be provided with the surface 
susceptibility in the MC results shown below. 




0.2 0.4 0.6 0.8 1 1.2 1.4 j 1.6 



FIG. 3: First two layer-magnetizations obtained by the Green 
function technique vs. T for J a = —0.5 with I = — I s = 0.1. 
The surface-layer magnetization (lower curve) is much smaller 
than the second-layer one. See text for comments. 

Figure 2] shows the non frustrated case where J s = 0.5, 
with I = I s = 0.1. As seen, the first-layer magnetization 
is smaller than the second-layer one. There is only one 
transition temperature. Note the difficulty for numerical 
convergency when the magnetizations come close to zero. 




FIG. 4: First two layer-magnetizations obtained by the Green 
function technique vs. T for J 3 = 0.5 with I = I s = 0.1. 

We show in Fig. [5] the phase diagram in the space 
(J S ,T). Phase I denotes the ordered phase with surface 
non collincar spin configuration, phase II indicates the 
collinear ordered state, and phase III is the paramagnetic 
phase. Note that the surface transition does not exist for 

Js > Jl 




FIG. 5: Phase diagram in the space (J a ,T) for the quantum 
Heisenberg model with N z = 4, / = \I S \ = 0.1. See text for 
the description of phases I to III. 



IV. MONTE CARLO RESULTS 

It is known that methods for quantum spins, such as 
spin-wave theory or Green function method presented 
above, suffer at high temperatures. Spin-wave theory, 
even with magnon-magnon interactions taken into ac- 
count, cannot go to temperatures close to T c . Green 
function method on the other hand can go up to T c but 
due to the decoupling scheme, it cannot give correct crit- 
ical behavior at T c . Fortunately, we know that quan- 
tum spin systems behave as their classical counterparts 
at high T . So, to see if the phase diagram obtained in the 
previous section for the quantum model is precise or not, 
we can consider its classical version and use MC simula- 
tions to obtain the phase diagram for comparison. MC 
simulations are excellent means to overcome approxima- 
tions used in analytic calculations for the high T region 
as discussed above. 

In this paragraph, we show the results obtained by MC 
simulations with the Hamiltonian ([1]) but the spins are 
the classical Heisenberg model of magnitude 5=1. 

The film sizes are Lx Lx N z where N z = 4 is the num- 
ber of layers (film thickness) taken as in the quantum case 
presented above. We use here L = 24, 36, 48, 60 to study 
finite-size effects as shown below. Periodic boundary con- 
ditions are used in the XY planes. The equilibrating time 
is about 10 6 MC steps per spin and the averaging time 
is 2 x 10 6 MC steps per spin. J = 1 is taken as unit of 
energy in the following. 

Let us show in Fig. [6] the layer magnetization of the 
first two layers as a function of T , in the case J s = 0.5 
with N z = 4 (the third and fourth layers are symmetric). 
Though we observe a smaller magnetization for the sur- 
face layer, there is clearly no surface transition just as in 
the quantum case. 

In Fig. [7] we show a frustrated case where J s = —0.5. 
The surface layer in this case becomes disordered at a 
temperature much lower than that for the second layer. 
Note that the surface magnetization is not saturated to 
1 at T = 0. This is because the surface spins make an 
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FIG. 6: Magnetizations of layer 1 (circles) and layer 2 (dia- 
monds) versus temperature T in unit of J/fcs for J s = 0.5 
with I — I 3 = 0.1, L = 36. 



FIG. 8: Susceptibilities of layer 1 (circles) and layer 2 (dia- 
monds) versus temperature T in unit of J/fcs for J s = 0.5 
with I = I S = 0.1, L = 36. 
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FIG. 7: Magnetizations of layer 1 (circles) and layer 2 (dia- 
monds) versus temperature T in unit of J/fcs for J 3 = —0.5 
with J = -J s = 0.1, L = 36. 



angle with the z axis so their z component is less than 1 
in the GS. 

Figures H] shows the susceptibilities of the first and sec- 
ond layers in the case where J s — 0.5 with I = I s = 0.1 
where one observes the peaks at the same temperature 
indicating a single transition in contrast to the frustrated 
case shown in Fig. [5] These results confirm the above 
results of layer magnetizations. 

To establish the phase diagram, the transition temper- 
atures are taken at the change of curvature of the layer 
magnetizations, i.e. at the maxima of layer susceptibil- 
ities shown before. Figure [10] shows the phase diagram 
obtained in the space (J S ,T). Interesting enough, this 
phase diagram resembles remarkably to that obtained for 
the quantum counterpart model shown in Fig. [5] 

Let us study the finite size effect of the phase transi- 
tions shown in Fig. 1101 To this end we use the histogram 
technique which has been proved so far to be excellent for 
the calculation of critical exponents ! 24 ' 25 ! 26 The principle 
is as follows. Using the Metropolis algorithm to deter- 
mine approximately the critical temperature region, then 
choosing a temperature To as close as possible to the pre- 
supposed transition temperature. We then make a very 
long run at Tq to establish an energy histogram. From 



FIG. 9: Susceptibility of layer 1 (circles) and layer 2 (dia- 
monds) versus temperature T in unit of J/ks for J s — —0.5 
with I = ~I S = 0.1, L — 36. Note that for clarity, the sus- 
ceptibility of the layer 2 has been multiplied by a factor 5. 
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FIG. 10: Phase diagram in the space (J S ,T) for the classical 
Heisenberg model with N z = 4, I = \I S \ = 0.1. Phases I to 
III have the same meanings as those in Fig. [5] . 



formulae established using the statistical canonical distri- 
bution, we can calculate physical quantities in a contin- 
uous manner for temperatures around Tb i 24 ' 25 ' 26 We do 
not have problem to identify the transition temperature 
as well as the maximal values of fluctuation quantities 
such as specific heat and susceptibility. 
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Figure [TT] shows the susceptibility versus T for L = 
36, 48, 60 in the case of J s = 0.5. For presentation conve- 
nience, the size L — 24 has been removed since the peak 
for this case is rather fiat in the scale of the figure. How- 
ever, it shall be used to calculate the critical exponent 
7 for the transition. As seen, the maximum ^ max of the 
susceptibility increases with increasing L. 
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FIG. 13: Layer susceptibilities versus T for L — 36, 48, 60 with 
J s — —0.5 and I = — I s — 0.1. Left (right) figure corresponds 
to the first (second) layer susceptibility. 



FIG. 11: Susceptibility versus T for L = 36, 48, 60 with J s 
0.5 and I = I s = 0.1. 
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FIG. 14: Transition temperature versus L for the surface layer 
in the case J a = —0.5 with I = — I s — 0.1. The inset shows 
the enlarged scale. 



FIG. 12: Transition temperature versus L with J a 
I = I s — 0.1. The inset shows the enlarged scale. 



0.5 and 



For completeness, we show in Fig. [12] the transition 
temperature as a function of L. A rough extrapolation 
to infinity gives T c °° ~ 1.86 ± 0.02. 

In the frustrated case, i.e. J s < Jf, we perform the 
same calculation for finite-size effect. Note that in this 
case there are two phase transitions. We show in Fig. [T3] 
the layer susceptibilities as functions of T for different L. 
As seen, both surface and second-layer transitions have 
a strong size dependence. 

We show in Fig. [H] and [15] the size dependence of 
the transition temperatures T\ (surface transition) and 
T2 (second-layer transition). 

The size dependence of the maxima observed above al- 
lows us to estimate the ratio j/f. We show now ln% max 
as a function of InL for the different cases studied above. 
Figure HHa) and (b) correspond respectively to the tran- 
sitions of surface and second layer occurring in the frus- 
trated case with J s = —0.5, while Figure fTBTc) corre- 
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FIG. 15: Transition temperature versus L for the second layer 
in the case J s = —0.5 with / = — I s — 0.1. The inset shows 
the enlarged scale. 



sponds to the unique transition occurring in the non frus- 
trated case with J s = 0.5. The slopes of these straight 
lines give 7 /V ~ 1.864 ± 0.034 (a), 1.878 ± 0.027 (b), 
1.801 ± 0.027 (c). The errors were estimated from the 
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FIG. 16: Maximum of surface-layer susceptibility versus L 
for L = 24,36,48,60 with J s = -0.5 (a,b), J s = 0.5 (c) 
and I — \I a \ — 0.1, in the In — In scale. The slope gives j/u 
indicated in the figure for each case. See text for comments. 



mean least-square fitting and errors on the peak val- 
ues obtained with different values of T (multi histogram 
technique). Within errors, the first two values, which 
correspond to the frustrated case, can be considered as 
identical, while the last one corresponding to the non 
frustrated case is different. We will return to this point 
later. 

At this stage, we would like to emphasize the following 
points. First, we observe that these values of "f/v are 
found to be between those of the two-dimensional (2D) 
Ising model (7/V = 1.75) and the three-dimensional one 
{l/ v — 1-241/0.63 ~ 1.97). A question which naturally 
arises on the role of the Ising-like anisotropy term, of the 
two-fold chiral symmetry and of the film thickness. The 
role of the anisotropy term and the chiral symmetry is 
obvious: the Ising character should be observed in the 
result (we return to this point below). It is however not 
clear for the effect of the thickness. Some arguments, 
such as those from renormalization group, say that the 
correlation length in the direction perpendicular to the 
film is finite, hence it is irrelevant to the criticality, the 
2D character therefore should be theoretically preserved. 
We think that such arguments are not always true be- 
cause it is difficult to conceive that when the film thick- 
ness becomes larger and larger the 2D universality should 
remain. Instead, we believe that that the finite thick- 
ness of the film affects the 2D universality in one way 
or another, giving rise to "effective" critical exponents 
with values deviated from the 2D ones. The larger the 
thickness is, the stronger this deviation becomes. The 
observed values of 7/1/ shown above may contain an ef- 
fect of a 2D-3D cross-over. At this point, we would like 
to emphasize that, in the case of simple surface condi- 
tions, i.e. no significant deviation of the surface param- 
eters with respect to those of the bulk, the bulk behav- 
ior is observed when the thickness becomes larger than 
a dozen of atomic layers^ surface effects are insignifi- 
cant on thermodynamic properties of the film. There are 
therefore reasons to believe that there should be a cross- 
over from 2D to 3D at some film thickness. Of course, 



this is an important issue which needs to be theoreti- 
cally clarified in the future. We return now to the effect 
of Ising anisotropy and chiral symmetry. The deviation 
from the 2D values may result in part from a complex 
coupling between the Ising symmetry, due to anisotropy 
and chirality, and the continuous nature of the classi- 
cal Heisenberg spins studied here. This deviation may 
be important if the anisotropy constant / is small. For 
the effect chiral symmetry, it is a complex matter. To 
show the complexity in determining the critical univer- 
sality with chiral symmetry, let us discuss about a simpler 
case with similar chiral symmetry: the XY model on the 
fully frustrated Villain's lattice. There has been a large 
number of investigations on the nature of the transition 
observed in this 2D case in the context of the frustration 
effects i 1 ' 27 ' 28 In this model, though the chirality symme- 
try argument says that the transition should be of 2D 
Ising universality class, many investigators found a de- 
viation of critical exponents from those of the 2D case 
(see review in Refi 2 ^). For example, the following values 
are found for the critical exponent v. v — U.889 2 -^ and 
v = 0.813^ These values are close to that obtained for 
the single transition in a mixed Al"-Ising model which is 
0.85 i 31 i 32 It is now believed that the XY character of the 
spins affects the Ising chiral symmetry giving rise to those 
deviated critical exponents. Similarly, in the case of thin 
film studied here, we do not deal with the discrete Ising 
model but rather an Ising-like Heisenberg model. The 
Ising character due to chiral symmetry of the transition 
at the surface is believed to be perturbed by the continu- 
ous nature of Heisenberg spins. The transition of interior 
layer shown in Fig. I16f b) suffers similar but not the same 
effects because of the absence of chiral symmetry on this 
layer. So the value 7/z/ is a little different. In the non 
frustrated case shown in Fig. [ToTc). the deviation from 
the 2D Ising universality class is less important because 
of the absence of the chiral symmetry. This small de- 
viation is believed to stem mainly from the continuous 
nature of Heisenberg spins. 

To conclude this paragraph, we believe, from physical 
arguments given above, that the deviation from 2D Ising 
universality class of the transitions observed here is due 
to, in an decreasing order of importance, the effect of the 
coupling between the continuous degree of freedom of 
Heisenberg spin to the chiral symmetry, the small Ising- 
like anisotropy and the film thickness. 



V. CONCLUDING REMARKS 

We have studied, by means of a Green function method 
and MC simulations, the Heisenberg spin model with an 
Ising-like interaction anisotropy in thin films of stacked 
triangular lattices. The two surfaces of the film arc frus- 
trated. We found that surface spin configuration is non 
collinear when surface antiferromagnetic interaction is 
smaller than a critical value J%. In the non collinear 
regime, the surface layer is disordered at a temperature 
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lower than that for interior layers ("soft" surface). This 
can explain the so-called " magnetically dead surface" ob- 
served in some materials . 33 i 34 The surface transition dis- 
appears for J s larger than the critical value J%. A phase 
diagram is established in the space (T, J s ). A good agree- 
ment between the Green function method and the MC 
simulation is observed. This is due to the fact that at 
high temperatures where the transition takes place, the 
quantum nature of spins used in the GF is lost so that 
we should find results of classical spins used in MC sim- 
ulations. We have also studied by MC histogram tech- 
nique the critical behavior of the phase transition using 
the finite-size effects. The result of the ratio of critical 
exponents j/u shows that the nature of the transition is 
complicated due to the influence of several physical mech- 



anisms. The symmetry of the ground state alone cannot 
explain such a result. We have outlined a number of 
the most relevant mechanisms. Finally, we note that in 
surface magnetism the low surface magnetization exper- 
imentally observe d 33 ! 34 has been generally attributed to 
the effects of the reduction of magnetic moments of sur- 
face atoms and/or the surface-localized low- lying magnon 
modes. The model considered in this paper adds another 
origin for the low surface magnetization: surface frustra- 
tion. It completes the list of possible explanations for 
experimental observations in thin films. 
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